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A)  Objectives 

Objectives  of  original  proposal  were  met  and  the  results  are  outlined  in  this  final  report. 

B)  Status  of  Effort  for  Duration  of  Funding 

In  the  funding  period,  several  breakthroughs  were  made  in  the  theoretical  study  and 
modeling  and  simulation  of  ionic  liquids  (ILs).  A  polarizable  model  was  developed  to  simulate 
ILs  more  accurately  at  the  atomistic  level.  A  multiscale  coarse-graining  methodology  was 
applied  to  IL  systems,  extending  computer  simulations  of  ILs  from  microscopic  to  mesoscopic 
scales.  With  these  advanced  simulation  methods,  a  spatial  heterogeneity  due  to  tail  aggregation 
was  discovered.  The  suggested  mechanism  explains  the  experimental  results,  and  changed  the 
way  people  think  of  IL  systems.  Based  on  the  tail  aggregation  phenomenon,  the  surface 
layering  of  IL  surface,  the  micelle  formation  in  IL/water  mixture,  and  the  aggregation  of 
hydrophobic  molecules  in  IL  systems  were  observed  and  systematically  explained.  They  have 
great  values  in  guiding  the  real  applications  of  ILs  in  various  conditions.  With  the  above  general 
understanding  of  IL  systems,  the  physical  properties  of  an  IL  candidate  for  propellant,  1- 
hydroxyethyl-4-amino-1,  2,  4-triazolium  nitrate  (HEATN),  were  studied  with  the  all-atom 
polarizable  model.  The  mechanism  suggested  for  HEATN  has  provided  deep  theoretical  insight 
for  this  specific  IL,  and  will  guide  more  efficient  storage  and  usage  of  HEATN  as  propellant. 


C)  Significant  Work  Accomplished 
1)  Atomistic  Polarizable  Model  Development 

Room  temperature  molten  salts  (RTMS),  or  ionic  liquids  (ILs),  stimulated  a  significant 
body  of  research.  Like  inorganic  molten  salts,  ILs  are  composed  solely  of  ions,  but  in  contrast 
their  melting  point  is  often  less  than  100  C.'  Despite  the  level  of  research  activity,  many  of  the 
properties  of  these  interesting  liquids  remained  elucidated,2  This  new  class  of  liquid  promoted 
more  and  more  experimental  and  computational  studies.  A  rich  level  of  detail  on  the  structural 
and  dynamical  properties  was  explored. 

For  an  IL,  it  was  found  both  in  simulation3'6  and  experiment,7  8  that  the  local  environment 
around  the  IL  ions  was  highly  anisotropic.  For  such  a  system,  the  particles  were  not  well 
described  as  rigid  bodies,  because  their  electron  densities  were  distorted  during  the  interactions 
with  each  other  and  the  changes  in  internal  structure  were  basically  a  many-body  effects,9  which 
made  it  highly  desirable  to  model  those  interactions  using  a  fully  polarizable  forcefield.'0  Given 
the  high  polarizability  of  IL  ions,11 12  it  was  a  key  priority  to  implement  a  polarizable  model  in  the 
computer  simulation  of  IL  systems. 

Importantly,  with  this  AFOSR-supported  project,  a  polarizable  forcefield  for  the  ionic 
liquids  such  as  l-ethyl-3-methylimidazolium  nitrate  (EMIMVNCV)  was  developed.12  The 
polarizable  model  was  developed  based  on  Thole's  smearing  dipole  model,13  in  which  each 
atomic  site  was  treated  as  a  polarizability  center  and  the  anisotropic  molecular  polarizability 
could  be  reproduced  through  the  interaction  among  the  atomic  sites. 13,14  Though  the  polarizable 
forcefield  was  initially  developed  for  the  IL  system  EMIMVNCV,  the  method  is  quite  general  and 
readily  applied  on  other  ILs,  such  as  the  1-hydroxyethyl-4-amino-1,  2,  4-triazolium  nitrate 
(HEATN). 


2)  Polarizable  Molecular  Dynamics  Code  Development 

In  conjunction  with  the  polarizable  IL  model  mentioned  above,  a  molecular  dynamics 
(MD)  Fortran  code  was  also  developed  and  implemented  in  the  DL_POLY  MD  simulation 
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package.15  The  long-range  muitipofe  interactions  (charge-charge,  charge-dipole,  and  dipole- 
dipole)  of  the  above  polarizable  forcefield  of  an  infinite  system  were  handled  with  the  Ewald 
summation  method.16,17  This  MD  code  for  the  polarizable  forcefield  was  also  reasonably  well- 
parallelized  with  the  replicated  data  method.16 

In  the  MD  simulations  with  the  polarizable  forcefield,  the  induced  dipole  could  be 
calculated  iteratively, 17,19  or  propagated  by  the  extended  Lagrangian  (ext-L)  method.20  Both 
methods  for  the  induced  dipole  were  implemented  in  the  MD  referenced  here.  For  the  iteration 
method,  the  induced  dipoles  are  solved  iteratively  and  generally  converged  in  5  -  6  cycles  with 
a  convergence  criterion  10'5  Debye.  For  a  system  size  of  -10,000  atoms,  the  iteration  method 
was  very  time  consuming,  even  when  run  on  multiple  processors.  With  the  ext-L  method, 
however,  the  induced  dipole  degrees  of  freedom  were  propagated  in  time,  with  a  proper  chosen 
fictitious  mass  for  each  dipole,  instead  of  calculating  them  iteratively.  Therefore,  just  one  cycle 


was  needed  for  the  ext-L  method. 

The  first  study  for  a  400  K  bulk 
simulation12  for  EMIM+/N03'  showed  that  the 
main  effect  of  the  electronic  polarizability 
would  make  the  ions  considerably  more  mobile 
when  compared  to  a  non-polarizable  model,  in 
which  the  electronic  polarizability  was  not 
explicitly  represented  in  the  forcefield.  This 
effect  resulted  in  a  much  lower  viscosity  and 
significantly  higher  diffusion  constant.  Fig,  1 
shows  the  mean  square  displacement  (MSD), 
on  a  log-log  plot  up  to  1  ns,  of  EMIM*  and  N03‘ 
for  both  the  polarizable  and  the  non- 
polarizable  models.  Remarkably,  the  slope 
of  the  MSD  of  the  polarizable  model  was  more 
than  three  times  larger  than  the  non- 
polarizable  model  at  the  end  of  1  ns.  The  short 
time  inertial  motion  due  to  ballistic  collisions 
lasts  less  than  0,1  ps,  followed  by  a  long  time 
intermediate  region  (-100  ps  for  the 
polarizable  model  and  -800  ps  for  the  non- 
polarizable  model),  finally  limited  to  linear 
diffusion.  A  linear  fit  to  the  MSD  between  900 
-  1000  ps  and  use  of  the  Einstein  relation, 
<Ar(t)z>=6Dt,  gave  self-diffusion  constants  of 
1 .49x1  O’10  m2/s  and  1.55x10‘10  m2/s  for 
polarizable  EMIM*  and  N03',  respectively,  and 
5.07x10'n  m2/s  and  4.80x1 0'11  m2/s  for  non- 
polarizable  EMIM*  and  N03‘,  respectively. 


Figure  I  Mean  square  displacements  of  EMIM+  and 
NOS- for  both  the  polarizable  and  non-polarizable 
models. 


3)  Multiscale  Coarse-Graining  Model  Development 

A  multi-scale  coarse-graining  (MS-CG)  approach21  was  applied  to  the  EMIM4/N03'  IL  to 
construct  a  MS-CG  model  for  ILs,  The  short-range  non-bonded  forces  between  CG  sites  were 
force-matched  by  an  extended  least-square  force-matching  method,  while  the  bonded  forces 
between  intramolecular  CG  sites  were  retrieved  from  the  configurational  data  from  the 
polarizable  atomistic  MD  simulations.  This  MS-CG  model  rebuilt  satisfactory  structural  and 
thermodynamic  properties  and  showed  significant  time  integration  efficiency  (factors  of 
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hundreds  speed-up).  Simulations  of  iarge  IL  systems,  for  long  periods  of  time,  are  therefore 
possible  with  this  MS-CG  model. 


4)  Spatial  Heterogeneity  and  Domain  Diffusion  in  Ionic  Liquids 


By  extending  the  MS-CG  model,  a 
spatial  heterogeneity  due  to  cationic  tail 
aggregation  has  been  discovered.22  It  has 
been  found  that,  with  sufficiently  long  side 
chains,  the  tail  groups  of  cations  aggregate  to 
form  spatially  heterogeneous  domains  {white 
spheres  in  Fig.  2).  This  is  understood  as  the 
result  of  competition  between  the  charged 
electrostatic  interactions  between  head  groups 
and  anions  and  the  collective  short-range 
interactions  between  the  neutral  tail  groups. 
This  aggregation  helps  to  explain  a  number  of 
experimentally23,211  observed  physical 
phenomena  in  ILs,  especially  the  mesoscopic 
liquid  crystal-tike  inhomogeneity,  which  exists 
in  IL  systems.25'28 

The  extended  MS-CG  model  was  then 
used  to  qualitatively  investigate  the  phase 
transition  of  tail  domains,29  It  was  found  that, 
with  a  sufficiently  long  alkyl  chain  at  a  low 
enough  temperature,  the  tail  domains  remain 
relatively  stable,  despite  the  diffusion  of 
individual  ions  in  the  liquid  phase;  With 
increasing  temperature,  the  average  tail 
domains  began  to  diffuse,  while  beyond  a 
transition  temperature,  their  average  density 
had  an  almost  uniform  distribution,  although 
the  tail  groups  continued  to  form  instantaneous 
domains.  Fig.  3  shows  the  average  tail  density 
of  the  BMIMVNO3'  system  with  64  ion  pairs 
calculated  by  the  newly  defined  heterogeneity 


Figure  2  Spatial  heterogeneity  due  to  tail  aggregation 
in  ionic  liquids,  (a)  All  CG  sites,  (b)  Tail  groups  only, 
(c)  Head  groups  only,  (d)  Head  groups,  tail  groups , 
and  anions. 


(a)  T=700K 


(b)  T=I100K 


o 


zj*  ft 

■'ii.o 

too 

6.0 

4  a 

'20 

'o.o 


Figure  3  Tail  domain  diffusion  in  ionic  liquids. 


order  parameter.  Each  sphere  represents 
the  ensemble-averaged  local  density  in  that 
location.  The  warmer  the  color,  the  higher  the 
tail  density.  It  can  be  seen  that,  below  the 
transition  temperature  Tc  =  1200  K,  the  tail 
domains  have  relative  stationary  positions; 
while  above  Tc,  they  move  around  so  that  on 
average  they  distribute  almost  uniformly. 

The  atomistic  polarizable  model  was  extended  to  the  XMIM+/N03'  ILs  with  the  carbon 
groups  on  their  cationic  side  chains  ranging  from  2  to  12.30  These  polarizable  atomistic  models 
were  used  to  refine  the  spatial  heterogeneity  discovered  earlier22  by  the  MS-CG  model.  While 
the  tail  groups  were  still  found  to  aggregate  (white  spheres  in  Fig.  2),  the  more  detailed 
atomistic  models  discovered  the  continuous  polar  network  composed  of  the  charged  cationic 
head  groups  and  the  anions  (yellow  and  red  spheres  in  Fig.  2),  as  reported  earlier  by  Lopes  and 
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Padua. 31  We  have  interpreted  the  formation  of  the  continuous  polar  network  as  the  following. 
Due  to  their  strong  electrostatic  interactions,  the  anions  always  stay  dose  to  the  head  groups. 
The  charged  groups  retain  their  local  structures  relatively  unchanged,  forming  a  continuous 
polar  network,  and  leading  to  different  global  structures  when  varying  the  side-chain  length.  The 
atomistic  simulation  results  are  quantitatively  comparable  to  experimental  results.  The  tail 
aggregation  mechanism  we  have  suggested  based  on  the  simulation  results  explains  very  well 
the  experimentally21,3032  observed  increasing  viscosity  and  decreasing  mass  density  with 
increasing  side-chain  length, 


5)  Roughness  and  Layering  of  Ionic  Liquid  Surface 


Figure  4  One  snapshot  of  ihe  ionic 
liquid  surface. 


The  IL/vacuum  interface  has  been  investigated  with 
both  the  atomistic33  and  the  MS-CG34  models.  The 
EMIM7N03'  IL  was  first  simulated  by  both  the  non- 
polarizable  and  the  polarizable  atomistic  models  at  T  =  400 
K.  The  cations  and  the  anions  were  found  to  have  an 
enhanced  density  on  the  interface.  By  using  the  second 
Legendre  polynomial,  it  was  clearly  shown  that  the  cations 
on  the  interface  tend  to  be  perpendicular  to  the  surface  with 
their  ethyl  tail  groups  protruding  outward.  One  snapshot  of 
the  IL  surface  is  given  in  Fig.  4.  Balls  and  sticks  represent 
the  atoms  on  the  surface,  while  lines  represent  the 
underlying  ones.  This  picture  clearly  shows  the  roughness 
of  the  IL  surface.  This  study  also  demonstrated  that  the 
polarization  effect  is  important  for  the  IL  surface.  In 
particular,  the  surface  tension  given  by  the  polarizable 
model  was  much  smaller  than  that  given  by  the  non- 
polarizabte  model  and  closer  to  that  of  the  experimental. 

The  simulation  results  by  the  MS-CG  model34  for  EMIM*/N03'  also  showed  the  same 
surface  behavior  of  cations.  In  order  to  study  the  effect  of  the  alkyl  side-chain  length  to  the 
surface  behavior,  the  side  chain  of  the  MS-CG  model  for  the  EMIM+/N03-  IL  was  elongated  to 
build  the  CG  models  for  XMIM+/N03-,  with  the  carbon  groups  on  their  cationic  side  chains 
ranging  from  4  to  12,  which  were  then  used  to  perform  the  surface  MD  simulations  for  the 
longer-chain  systems.  The  surface  tension 
decreased,  while  the  side-chain  length 
increased.  A  multilayer  ordering  appeared 
more  clearly  for  longer  cationic  side-chain 
systems  (shown  in  Fig.  5).  Detailed 
analysis  indicated  that  such  multilayer 
ordering  is  also  resulted  from  the  strong 
short-range  interaction  of  the  tail  groups. 

Thus,  the  behavior  of  the  IL  surfaces  were 
governed  by  the  tail  aggregation 
mechanism  in  a  two  dimensional  regime. 

This  mechanism  better  explains  the 
experimental  results35"39  than  does  the 
mechanisms37,3*  suggested  by  others. 

Expected,  is  the  unification  of  the 

understanding  of  the  surface  behavior  of  IL  „.  ,  ,,  .  ,  , 

systems  and  further  Bmded  experimental  Jmalltlrtc„mPallLsur^“ 
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and  theoretical  investigations  of  IL  interfaces. 


6}  Micelle  Formation  in  Ionic  Liquid/Water  Mixture 


MD  simulations'*0  with  the  non-polarizable 
atomistic  force  field  were  carried  out  to  investigate 
the  influence  of  water  on  the  spatial 
heterogeneous  distribution  of  ionic  liquids.  It  was 
found  that,  by  increasing  the  amount  of  water,  the 
spatial  heterogeneity  of  the  ionic  liquid  first 
increases,  and  then  stays  on  its  maximum  in  a 
narrow  range  of  water  concentration,  before  it 
decreases  drastically.  Fig.  6  shows  one  snapshot 
for  the  IL/water  mixture  with  a  water  molar  fraction 
of  0.75,  when  the  IL  exhibits  the  strongest 
aggregation.  The  cationic  tail  groups  (yellow 
spheres)  aggregate  to  form  relatively  isolated  tail 
domains,  bordered  by  the  continuous  polar 
network  composed  of  the  cationic  head  groups 
(red  spheres)  and  the  anions  (pink  spheres),  while 

the  water  molecules  (blue  spheres)  are  embedded  in  the  polar  network.  In  contrast,  the  ionic 
diffusion  monotonically  grows  with  the  aqueous  concentration.  These  phenomena  are 
understood  as  the  result  of  the  continuous  competition  between  the  enhanced  hydrophobic 
interaction  of  the  cationic  tail  groups  and  the  attenuation  of  the  ion-ion  interactions  with 
increasing  amount  of  water,  This  mechanism  successfully  explains  the  experimentally'*1'42 
observed  micelle  formation  of  ionic  liquids  in  aqueous  solutions  and  agrees  with  the  general 
surfactant  critical  micelle  concentration  theory.43 


Figure  6  Micelle  Formation  in  !U water  mixture. 


7)  Aggregation  of  Hydrophobic  Molecules  in  Ionic  Liquids 


Continued  MD  simulations  were 
carried  out  to  investigate  the  aggregation 
of  small  non-polar  particles  in  ionic 
liquids.44  Fig.  7  shows  the  aggregation  of 
methane  molecules  in  BMIM+/CL‘.  The 
aggregation  of  non-polar  molecules 
decreases  with  increasing  symmetries  of 
the  anions,  or  equivalently,  decreases  with 
decreasing  strength  of  electrostatic 
attraction  between  polar  head/anion 
groups.  For  a  given  anion,  the  aggregation 
decreases  with  increasing  length  of  the 
side  chain.  A  comparison  between 


Figure  7  Aggregation  of  hydrophobic  molecules  in  ionic 
liquids. 


methane  and  butane  shows  that  methane  molecules  have  much  stronger  aggregation  tendency 
than  butane.  The  above  results  can  be  explained  by  a  unified  entropic  origin  within  the 
framework  of  the  previous  revealed  spatial  heterogeneity  phenomenon  due  to  tail 
aggregation22-29,30. 
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8)  Physical  Properties  of  HEATN 


An  atomistic  polarizable  model  was  developed  for  l-hydroxyethyl-4-amino-l,  2,  4- 
triazolium  nitrate  (HEATN),  an  energetic  IL,  which  is  a  good  candidate  for  ionic  fuel  and 
propellant,  and  used  in  the  MD  simulations  to  investigate  its  physical  properties.  The  upper 
panels  in  Fig,  8  show  the 


mean  square  displacements 
of  ions  with  different  time 
scales  at  various 
temperatures.  The  two  lower 
panels  in  Fig.  8  show  the  self¬ 
diffusion  coefficients  vs. 
temperature,  A  glass 
transition  point  at  T  -  275  K 
was  unambiguously  identified. 

A  complex  hydrogen  bond 
network  was  revealed  in 
HEATN,  which  is  responsible 
for  the  poor  packing  efficiency 
of  ions  and  the  low 
melting/glass  transition  point 
of  HEATN.  The  densities 
calculated  from  the  simulation 
data  are  1 .436  g/cm3  and 
1.427  g/cm3,  respectively,  for 
the  polarizable  and  non- 

polarizable  models.  The  viscosities  were  fitted  as  64.25  and  38.21  c.p.,  for  the  non-polarizable 
and  polarizable  models,  respectively.  Essentially,  both  the  high  viscosity  and  poor  ion  packing  of 
HEATN  can  be  attributed  to  the  introduction  of  hydroxyl  and  amino  groups  in  its  molecular 
structure. 


200  2S0  300  350  400  450 


T/K 

Figure  S  Mean  square  displacements  and  diffusion  constants  of  HEATN  at 
different  temperatures. 


9)  Extendable  Multiscale  Coarse-Graining  Models  Development 


The  efficient  use  of  ILs  required  the 
technique  of  systematic  design  of  ILs.  The 
analysis  of  the  interactions  between  building 
blocks  was  essential  for  the  development  of  this 
technique.  The  MS-CG  methods  developed  can 
obtain  the  effective  forces  between  CG  sites 
(building  blocks). 

An  effective  force  coarse-graining 
(EFCG)  method45  was  developed  as  a 
complement  to  the  general  MS-CG  method  and 
applied  to  the  EMIMVN03*  IL.  By  compromising 
some  accuracy  of  the  local  structures,  the  EFCG 
force  fields  were  transferable  between  different 
thermodynamic  conditions.  The  EFCG  force  field 
for  EMIMVNO3'  was  found  to  be  applicable  when 
generating  the  force  fields  for  similar  ILs  with 


Figure  9  Radial  Distribution  Functions  between  tail 
groups  for  the  IL  systems  with  different  alkyl  side- 
chain  lengths. 
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longer  alkyl  side  chains  by  rescaling  the  effective  electrostatic  forces  according  to  the  partial 
charges  on  the  CG  sites.46  Fig.  9  demonstrates  the  radial  distribution  functions  between  the 
cationic  tail  groups  obtained  from  the  CG  MD  ran  with  the  extended  EFCG  force  fields  for  the 
iLs  with  different  alkyl  side-chain  lengths.  The  numbers  represent  the  number  of  carbon  groups 
on  the  side  chain.  Fig.  9  clearly  shows  that  when  increasing  alkyl  side-chain  length,  the 
aggregation  of  the  cationic  tail  groups  became  stronger,  and  this  tendency  agreed  with  that 
observed  in  previous  all-atom  MD  simulations.30  Because  the  EFCG  force  fields  were 
extendable  for  longer  side-chain  IL  systems  and  transferable  between  different  thermodynamic 
conditions,  they  were  especially  suitable  for  constructing  the  effective  force  library  of  IL  building 
blocks.  Choosing  building  blocks  of  ILs  according  to  the  features  of  their  effective  forces  allowed 
systematic  design  of  ILs  to  meet  customized  application  requirements. 
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C.  Executive  Summary: 

1)  Production  of  Thesis  Dissertations 

None 

2)  Significant  Advances  and  Conclusions 

The  focus  of  this  AFOSR  project  was  the  molecular  dynamics  (MD)  simulation  of  ionic 
liquid  (IL)  structure,  dynamics,  and  interfacial  properties,  as  well  as  multiscale  descriptions  of 
these  novel  liquids  (e.g,,  bridging  with  fluid  mechanics  methods).  The  overall  importance  of  ionic 
liquids,  as  well  the  challenges  for  the  future  in  this  area,  were  at  the  time  of  this  proposal  and 
continue  to  be,  dearly,  high  priorities  for  the  Air  Force  advanced  propulsion  program.  During  the 
course  of  award  period,  MD  simulations  were  employed  to  study  the  physical  properties  of  an 
energetic  ionic  liquid,  1-hydroxyethyl-4-amino-1,  2,  4-triazolium  nitrate  (HEATN),  The 
aggregation  of  small  non-polar  molecules  have  been  observed  to  strongly  aggregate  in  the  ILs, 
and  the  strength  of  aggregation  varies  for  different  cation/anion  combinations.  A  transferable 
and  extendable  multiscale  coarse-graining  (MS-CG)  model  was  developed  for  the  ILs,  which 
enabled  the  systematic  design  of  ILs  at  the  CG  level,  When  combining  accomplishments,  which 
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includes  developing  an  all-atom  polarizable  model  for  ILs,  developing  the  MS-CG  models  for 
lls,  discovering  the  spatial  heterogeneity  in  ILs  due  to  tail  aggregation  and  the  tail  domain 
diffusion  phenomena,  revealing  the  surface  layering  phenomenon  on  the  IL  surface,  and 
investigating  the  influence  of  water  concentration  on  the  tail  aggregation,  the  research  efforts 
lead  to  significant  new  predictions  concerning  the  fundamental  behavior  of  ILs  as  a  function  of 
their  chemical  composition  and  open  the  door  towards  systematic  design  of  ILs  to  meet 
customized  application  requirements, 
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Chinese  Theoretical  and  Computational  Chemistry  Conference,  Hong  Kong,  January, 
2005. 


3.  ‘‘Multi-Scale  Simulations  of  Ionic  Liquids,"  AFOSR  Contractor’s  Meeting,  Monterey 
California,  May,  2005. 

4.  "Systematic  Coarse-graining  of  Condensed  Phase  Systems,"  Materials  Division 
Seminar,  Air  Force  Research  Laboratory,  Wright-Patterson  Air  Force  Base.  December, 
2005. 
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Comprehensive  Final  Report:  Gregory  A.  Voth.  PI 

ConiracVGrant  title:  Multrseale  Simulation  of  High  Energy  Density  Ionic  Liquids 
Conlr act/Grant#:  FA9550-01-1-0082 

Reporting  Period:  Total  Award  Period:  1 1/18/2003  —  12/31/2006 


5.  “Molecular  Dynamics  and  Coarse-grained  Simulations  of  Ionic  Liquids,"  Invited  Seminar 
at  the  Air  Force  Office  of  Scientific  Research  Workshop  on  Ionic  Liquids,  University  of 
Alabama,  February,  2006. 

6.  “Multiscale  Simulation  of  Ionic  Liquids:  An  Emerging  Picture  of  Their  Unique  Dynamical 
and  Structural  Behavior,"  Invited  Seminar  in  the  Symposium  on  Ionic  Liquids:  Not  Just 
Solvents  Anymore,  ACS  Division  of  Industrial  and  Engineering  Chemistry,  National 
American  Chemical  Society  Meeting,  Atlanta,  Georgia,  March  2006. 


(b)  Consulting/Advisory  Functions: 

Gaussian,  Inc, 

(c)  Transitions: 

The  simulation  capability  developed  in  this  project  is  available  to  DoD  researchers  at 
Edwards  AFB  and  elsewhere  to  aid  in  the  design  and  analysis  of  new  ionic  liquids  for 
propulsion  purposes. 

New  Discoveries,  Inventions,  Patents: 

None 

Honors/Awards: 

John  Simon  Guggenheim  Memorial  Fellowship,  2004  -  2005. 
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